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ABSTRACT 

Since the discovery of kHz quasi-periodic oscillations (QPO) in neutron star 
binaries, the difference between peak frequencies of two modes in the upper part 
of the spectrum, i.e. Au = uj h — uj k has been studied extensively. The idea that 
Auj is constant and (as a beat frequency) is related to the rotational frequency of 
the neutron star has been tested previously. The observed decrease of Au; when 
Uh and Uk increase has weakened the beat frequency interpretation. We put 
forward a different paradigm: a Keplerian oscillator under the influence of the 
Coriolis force. For such an oscillator, and the assumed Keplerian frequency 
ojk hold an upper hybrid frequency relation: lo\ — uj\ = 4Q 2 , where Q is the 
rotational frequency of the star's magnetosphere near the equatorial plane. For 
three sources (Sco X-l, 4U 1608-52 and 4U 1702-429), we demonstrate that the 
solid body rotation Q = Q Q = const, is a good first order approximation. Within 
the second order approximation, the slow variation of Q as a function of ujk 
reveals the structure of the magnetospheric differential rotation. For Sco X-l, 
the QPO have frequencies ~ 45 and 90 Hz which we interpret as the 1st and 
2nd harmonics of the lower branch of the Keplerian oscillations for the rotator 
with Q, not aligned with the normal of the disk: ujl/2t: = (Q/n)(uK/u}h) sin S 
where 5 is the angle between f2 and the vector normal to the disk. 

Subject headings: accretion, accretion disks — radiation mechanisms: thermal 
— stars:neutron — X-rays: general 
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1. Introduction 

Kilohertz quasi-periodic oscillations (KQPO) have been discovered by the Rossi X-ray 
Timing Explorer (RXTE) in a number of low mass X-ray binaries (Strohmayer et al. 1996). 
The existence of two observed peaks with frequencies uj k and uj h in the upper part of the 
QPO spectrum became a natural starting point in modeling the phenomenon. Attempts 
have been made to relate ujk, 0Jh and the peak difference frequency Au = Uh — ujk with the 
neutron star spin and possible Keplerian motion of the hot matter surrounding the star [see 
origin of Keplerian motion discussed in Titarchuk, Lapidus & Muslimov 1998 (hereafter 
TLM) and Titarchuk & Osherovich 1999]. 

In the beat frequency model, the burst oscillations and the kHz peak separation Ao; 
are both considered to be close to the neutron star spin frequency and thus Aa; is predicted 
to be constant (see, e.g., review by van der Klis 1998). However, recent observations of kHz 
QPO's in Sco X-l, 4U 1608-52 and 4U 1702-429 showed that Aa; decreases systematically 
when uo K and oo h increases (van der Klis et al. 1997; Mendez et al. 1998; Markwardt, 
Strohmayer & Swank 1999). Indications of such a decrease was found by Wijnands et al. 
(1998) and Ford et al. (1998) for 4U 1735-44 as well. Psaltis et al. (1998) showed that the 
measured Ao; in all other known sources is consistent with the values found in Sco X-l and 
4U 1608-52. These sources have become standards for comparison between the kHz QPO 
phenomena in different sources because of the high accuracy of measurements and the wide 
range of values of uj k and Uh available. For Sco X-l and 4U 1608-52, Au/2n changes from 
320 Hz to 220 Hz when u K /2ir changes from 500 Hz to 850 Hz 

In the lower part of the Sco X-l spectra, van der Klis et al. (1997) found two frequencies 
~ 45 Hz (referred to below as ul/2tt) and 90 Hz which also slowly increase with increase 
of ujk and Uh', presumably 90 Hz is the second harmonic of ul/^k. Any consistent model 
faces a challenging task of describing the dependence of Aa; and Wi on ojk and u h for the 
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available sources. 

In the beat-frequency model (BFM) considered by Miller, Lamb & Psaltis (1998), the 
kHz QPO's are identified with the orbital frequency at the sonic radius and its beat with 
the neutron star spin, and the lower frequency of horizontal branch oscillations (HBO) in 
Z-sources with the beat of the orbital motion at the magnetospheric radius with the neutron 
star spin. In the general-relativistic (GR) precession/apsidal motion discussed by Stella 
and Vietri (1998,1999), the kHz peaks are due to orbital motion and GR apsidal motion 
of a slightly eccentric orbit, and low-frequency QPO to the Lense-Thirring precession of 
the same orbit. Titarchuk and Muslimov (1997) suggested that a set of peaks seen in the 
kHz QPOs can be naturally explained in terms of the effect of rotational splitting of the 
main oscillation frequency in the disk. The type of peak (eigenmodes) is dependent on 
mode numbers and a parameter describing disk structure. In TLM the authors argue that 
this effect of the rotational splitting originated in a disk region identified as the centrifugal 
barrier (CB) region oscilatted in the vertical and radial directions. The size oscillations of 
the CB disk region and reprocessing of the disk photons there (due to Comptonization) led 
to an oscillation signal. 

In this Letter, we take a different approach and we suggest that a kHz QPO can be 
explained as oscillations of large scale inhomogeneities (hot blobs) thrown into the neutron 
star's magnetosphere. Participating in the radial oscillations with Keplerian frequency u>k, 
such blobs are simultaneously under the influence of the Coriolis force. The present model 
is very different than the TLM one where the site of the rotational splitting was associated 
with the CB disk region. Thus, studying the Keplerian oscillator in the rotating frame of 
reference (magnetosphere), we attempt to relate Acj and ujl with ujk and u>h- In our model, 
we compare results of linear theory with the observations for three sources. We stress the 
fundamental role of the Keplerian frequency (which we identify with uj k ) and the upper 
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hybrid frequency relation between lu k and u^. In this first Letter, we depict the main 
features of the kHz QPO phenomena. The discussion identifies some limitations of our 
approach and suggests further directions of the research we shall pursue. 



2. Keplerian Oscillator Under the Influence of the Coriolis Force 

We assume that the interplay between the centrifugal, gravitational force and magnetic 
force maintains radial oscillations (along the axis x) with the Keplerian frequency 



u K 



where G is the gravitational constant, M is the mass of the compact object and R is the 
radius of an orbit. Self-similar sp hero mak- type models of a magnetic atmosphere can serve 
as an example of the exact time dependent MHD solutions in which non-radial equilibrium 
is maintained by magnetic force, while all motion is strictly radial (Farrugia et al. 1995 
and references there in). Taking the magnetospheric rotation with angular velocity fl (not 
perpendicular to the plane of the equatorial disk), the small amplitude oscillations of such a 
blob thrown into the magnetosphere are described by equations (Landau and Lifshitz, 1960) 

x + uj 2 k x - 2yVl cos 5 = (2) 

y + 2xQ cos 5 — 2zVL sin 5 = (3) 

z + 2yQ sin 5 = (4) 

where 5 is the angle between O and the vector normal to the plane of radial oscillations, 
y is the azimuthal component of the displacement vector and z is the vertical component 
which is perpendicular to the plane of Keplerian oscillations. The dispersion relation for 
the frequency uj 

uj 2 [lo 4 - (uj 2 k + 4fi 2 V 2 + 4fi 2 c4 sin 2 5} = (5) 
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besides the non-oscillating mode (u — 0), describes two spectral branches (eigenmodes), 
high and low: 



u ; 1,2 * y°) 



ui 2 h ±(ui-im 2 u 2 K sin 2 <J)V2 

^ J 1,2 - - 

where 



c^ = (c4 + 4fi 2 ) 1 / 2 (7) 

is the analog of the upper hybrid frequency f uh in plasma physics (Akhiezer et al., 1975; 
Benson 1977). When the angle 8 is small (as we believe is indeed our case): 

iO X = LU h (8) 

is related with the radial eigenmode and 

uj 2 =0J L = 2Vt{ijj K /u} h )sin5 (9) 

is related with the vertical eigenmode. For 5cl, corrections for uji are of the second order 
in 5, but for u>2, the angle dependent term is the main term. 



3. Comparison of the Model with Observations 

According to formula (7) 

a-M-f 1 " (10 ) 

is the rotational frequency, which, for the highly conductive plasma in the case of corotation, 
should be approximately constant (Q = Qq). Interpreting the two frequencies in the upper 
part of the kHz spectrum as ujk and u>h [where ojk < 0Jh\ for Sco X-l, we plotted in 
Figure 1, v = Q/2n versus v K = uj k /2n (solid circles). Indeed, one may notice that as 
a first approximation Q = Q = const. The slow (but systematic) variation of Q/2it is 
between 330 and 350 Hz which should be compared with the variation by a factor of 1.5 for 
Au = Uh — ujk for the same range of ujk- If the magnetosphere corotates with the neutron 
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star (solid body rotation), then our procedure would determine the spin rotation of the 
star. In fact, there is a differential rotation profile which depends on the magnetic field B 
which can be presented in terms of the Chandrasekhar potential A (Chandrasekhar, 1956) 

p _ 1 ( IdA 1 OA B*(A) \ 

R\ Rdfji' (1 - /i 2 )V2 dR' (l-fi 2 ) 1 / 2 J { ' 

where \x = cos 6, 6 is the colatitude in the spherical system of coordinates and B* (A) is a 
function of A only. For a combination of dipole, quadrupole and octupole fields 



A = A (l-/i 2 ) 



5 + ^"^ a + 2 ) 



(12) 



where A , q and b are constants. For a highly conductive plasma, the MHD induction 
equation leads to the iso-rotation theorem of Ferraro (see Alfven and Falthammar, 1963) 
which in our notation means that Q is a function of A only. Following a recent suggestion 
of Osherovich and Gliner (1999), we assume that in the vicinity of a magnetic star Q(A) 
has a maximum and therefore can be presented as 

Q(A) = n - a 2 A 2 (13) 

where Qo and a are constants. Since in the equatorial plane (/i = 0), according to (12), the 
quadrupole term does not contribute to A, from formula (13) we find 

^-"•-(i-s)' (14) 

where a' and (3' are positive constants. Expressing R through the Keplerian frequency 
v K = uo k /(2t[) from equation (1), we find that formula (14) leads to the following 
presentation for Q: 

Q(u K )/2n = C + C x vf + C 2 uf + C 3 4 (15) 

where C 2 = 2\fC\C^ (not a free parameter because the radial part of the expansion is a full 
square of z/^ 3 and v\ with new positive constants a and (3). The fit for Sco X-l (solid line 
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in Figure 1) by formula (16) demonstrates the agreement of the proposed model with the 
data. Indeed, for Sco X-l, we find C = n /2n = 345 Hz, C x = -a 2 = -3.29 • 1(T 2 Hz" 1 / 3 , 
C 2 = 2a{3 = 1.017- 10" 5 Hz" 5 / 3 and C 3 = -[3 2 = -7.76- lO" 10 Hz" 3 with high X 2 = 37.6/39. 

In our fitting procedure, C , Ci, C 2 and C 3 are treated as four free parameters, i.e., 
the relation between C 2 and the two constants C\ and C 3 is not enforced. The fact that the 
derived C 2 indeed satisfies the relation C 2 = l^fC^Cl shows the internal consistency of the 
model. 

Three available data points for the 4U 1702-429 source are consistent with a constant 
value of rotation frequency Q /2n = 380 Hz, but the small number of data points does not 
allow us to recover the profile of Q(vk)- However, the fit of 4U 1608-52, shown in Figure 
2, demonstrates that the excellent agreement between our model and observations for Sco 
X-l is not accidental. For 4U 1608-52, we found C = tt /2ir = 345 Hz, d = -3.91 • 10~ 2 
Hz 1 / 3 , C 2 = 1.23 ■ 10~ 5 Hz" 5 / 3 and C 3 = -9.45 • 10" 10 Hz~ 3 with X 2 = 12.02/13. 

Sco X-l allows us to check the prediction of the model for oo L [formula (10)]. Taking the 
profile VL{i>k), presented in Figure 1, we plot u L = u L /2n in Figure 3. Choosing 5 = 5.5°, we 
found that our plot for v L and 2u L fit the data for the observed frequencies of 45 Hz and 
90 Hz, respectively. From these results, we conclude that the rotational axis of Sco X-l and 
the normal to the plane of radial oscillations do not coincide but are separated by 5.5°. 



4. Summary and Discussion 

The main thrust of this work is based on the hypothesis that two of the observed 
frequencies in the upper part of the kHz QPO spectrum are the Keplerian frequency 
(wjc) and the upper hybrid frequency (u>h) of the Keplerian oscillator in the frame of 
reference rotating with the angular frequency fl For three sources, we showed that 
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calculated according to the formula f2 = — u 2 K /2 changes with ujk significantly 
less than Auj = Uh — ujk, suggesting approximately constant rotation of the star's 
magnetosphere. The detailed profile of fl is successfully modeled as a function of u K within 
the dipole-quadrupole-octupole approximation of the magnetic field for Sco X-l and for 
4U 1608-52. For Sco X-l, the derived 0,(u K ) is used to compare predictions of the model 
for the lower branch frequency u>l/27t with observations of 45 Hz and 90 Hz (presumably 
2uj l /2tx). The last fit is done with only one unknown parameter 5 (angle between the fi 
and the normal vector to the plane of Keplerian oscillations). 

A good fit for uj l and 2u L confirmed the suggested model and revealed the expectedly 
small angle 8 = 5.5°. The profile for £l(v K ) has been modeled strictly in the equatorial 
plane (/i = 0) for the ideally aligned rotator (5 = 0). With 5^0, the next approximation 
will reveal the asymmetry of the QPO spectrum due to the quadrupole term. It has 
been shown (Osherovich, Tzur and Gliner, 1984) for the solar corona that the quadrupole 
term introduces North-South asymmetry in the makeup of the magnetic atmosphere 
surrounding a gravitating object. We shall pursue a similar study for the neutron star 
with the expectation of finding a signature of such asymmetry in the kHz QPO. Recent 
measurements show that the octupole field is comparable to the global dipole field near the 
sun, while the quadrupole field contributes 20 — 30% of the total field (Wang et al., 1997; 
Osherovich et al., 1999). 

In this Letter, we restrict the scope of the work to a model for the linear Keplerian 
oscillator. The presence of the second harmonic ( 90 Hz) of the 45 Hz mode, noticed by van 
der Klis et al. (1997 ), strongly suggests that a weakly nonlinear model is desirable. Some 
oscillations below uo L (for Sco X-l oscillations with frequencies 10 — 20 Hz) we attribute 
to the innermost edge of the Keplerian disk, which adjusts itself to the rotating central 
object (i.e., neutron star). The physics of these oscillations related to accreting matter with 
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a large angular momentum is outside the scope of this Letter and it have been described in 
a separate paper (Titarchuk & Osherovich 1999). 

We are grateful to J. Fainberg and R.F. Benson for help in the preparation of this 
paper and discussions and for fruitful suggestions by the referee. L.T. thanks NASA for 
support under grants NAS-5-32484 and RXTE Guest Observing Program and Jean Swank 
for support of this work. In particular, we are grateful to Michiel van der Klis and Mariano 
Mendez for the data which enabled us to make detailed comparisons of results of our model 
with the data. 
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Fig. 1. — The inferred rotational frequency of the NS magnetosphere (see Eq. 10) as a 
function of the frequency of the low kHz QPO peak (Keplerian frequency). For Sco X-l 
(van der Klis et al. 1997) (solid circles) and for 4U 1702-429 (Markwardt et al. 1999) (open 
squares). The solid (in the case of Sco X-l) line and dashed line (in the case of 4U 1702-429) 
are theoretical curves calculated using the multipole expansion of the rotational frequency 
(Eq.15 ). 



Fig. 2. — The inferred rotational frequency of the NS magnetosphere (see Eq. 10) as a 
function of the frequency of the low kHz QPO peak (Keplerian frequency). For for 4U 1608 
(Mendez et al. 1998) (solid squares). The solid line is the theoretical curve (see Eq.15). 



Fig. 3. — HBO frequency vs high peak KQPO frequency, for Sco X-l (van der Klis et al. 
1997) (solid circles). The solid line is theoretical dependence of the low branch oscillation 
frequencies (40-50 Hz) and second harmonics (80-100 Hz) on the high hybrid frequencies, 
calculated using the multipole expansion of the rotational frequency (Eq. 16) with 5 = 5.5°.. 



